function f=density_X(X,x_interest,kernel_type,h_x)

N=size(X,1);
K=size(x_interest,1);
term_x=kfunc((X*ones(1,K) - ones(N,1)*x_interest')/h_x,kernel_type);
f=(1/(N*h_x))*sum(term_x,1);

    %


